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Abstract 

In 1970 Baxter considered the statistical three-coloring lattice model for the case 
of toroidal boundary conditions. He used the Bethe ansatz and found the partition 
function of the model in the thermodynamic limit. We consider the same model 
but use other boundary conditions for which one can prove that the partition func- 
tion satisfies some functional equations similar to the functional equations satis- 
fied by the partition function of the six-vertex model for a special value of the 
crossing parameter. 



1 Introduction 



We consider various colorings of a square n x m lattice with three colors, such that 
any two adjacent faces have different colors. Sometimes one imposes some boundary 
conditions restricting possible colorings. The most common case here is the toroidal 
boundary conditions, when the colors of the first face must be different from the color 
of the last one for each row and column of the lattice. An example is given in Figure Q3 
where we use a convenient labeling of the colors by the elements of the ring Z3. We 



Figure 1: 

denote the elements of Z3 by 0, 1 and 2. When an element of Z3 arises in a context 
where an integer should be, it is treated as the corresponding integer 0, 1 or 2. Vice 
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verse, when an integer i arises in a context where an element of Z3 should be, it is 
treated as the element of Z3 corresponding to the reminder after division of i by 3. It is 
clear that the colors of the adjacent faces are different if and only if the corresponding 
elements of Z3 differ by +1 or —1. 

The simplest combinatorial problem arising here is to enumerate all possible col- 
orings. More generally one can try to find the numbers C ntm {k.Q,kykj) of colorings 
with k-Q faces of color 0, kj faces of color 1, and k^ faces of color 2. Certainly, the num- 
bers C nim (kQ, ky kj) are different from zero only if k^ + kj + kj = nm. It is useful to 
introduce the generating function 



where Zq, Zj, Zj are arbitrary numbers. 

One can treat colorings as states of a quantum statistical model regarding colors as 
states of faces. Here the numbers Zq, z v z^_ are the corresponding Boltzmann weights 
of faces, the Boltzmann weight of a state is the product of the Boltzmann weights 
of the faces, and Z„ /ffZ (z^, Zy, z^) is the partition function (state sum) of the model. It 
is convenient to represent the weight of a state as the product of the weights of the 
vertices of the lattice. Here if a vertex belongs to four faces of colors r\, r^, r^, and r^, 
we assign to it the weight (z ri z f - 2 z r3 z f - 4 ) 1/4 . One can generalize the model assuming 
that the weight of a vertex is an arbitrary function of the colors of the four adjacent 
faces. It is also possible to introduce a dependence on a spectral parameter and satisfy 
the Yang-Baxter equations (star-triangle relation) to obtain an integrable system. 

In 1970 Baxter, using an appropriate version of the Bethe ansatz, found the partition 
function Z M/?M (zq, Zy, z^) for the toroidal boundary conditions in the thermodynamic 
limit [HI3. After that, one of the authors of the present paper constructed a solution of 
the Yang-Baxter equations which, for the special value of the spectral parameter, leads 
to the partition function of the three-coloring model [3J. With the help of the inversion 
trick [4], the largest eigenvalues for the family of commuting transfer matrices were 
found and the results of the paper [1J were reproduced. 

To describe the solution of the Yang-Baxter equations given in the work [3], note 
that there are six types of possible vertex color configurations given by the first raw 
of the pictures in Figure |2l where r is one of the colors 0, 1, and 2. The corresponding 
weights of the work ||3]] are expressed via standard elliptic theta functions of nome 
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as follows 
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Figure 2: 



y{cp\rj) j'((p\n) 



Here cp is the spectral parameter associated with the central vertex of a four-face con- 
figuration, A is an arbitrary fixed parameter, and 



Cr(A,p) 



4 (A + 2n(r - l)/3|p)0 4 (A + 2n(r + l)/3\p) 



el(A + 2nr/3\p) 
It is useful to have in mind thafl 

and that 

a r+ i(<p|A) = a r (^|A + 27r/3), 
)6 r+T (<p|A) = a r (^|A + 27T/3), 
7 r+ T(<pl A ) = 7r(<p|A + 27T/3), 

For cp = we have 

« r (0) = <(0) = I) 1 *, 



= a^(^|A + 27T/3), 
# +T (?|A)=#(?|A + 27r/3), 
7j +T (?|A) =7^|A + 27r/3). 



0,(0) = #(0) = # /4 , 

7;(o) = ^ 2 T ^ i/2 , 



(5) 

(6) 

(7) 
(8) 
(9) 



and, using the equality ©, it is not difficult to see that we come to the Baxter's three- 
coloring model with 

It follows from the equality © that ZqZjZ^ = 1, but this fact does not lead to any loss 
of generality 

In 1961 Lenard remarked that there is a correspondence between the three-colo- 
rings and the states of the six-vertex model [6J. Recall that the six-vertex model is 
defined on a square lattice. A state of the model is specified by a choice of the direction 
of each internal edge by placing an arrow on it. The arrows obey the rule, called the 



1 Below, when it does not lead to misunderstanding, we do not write up explicitly dependence on 
the fixed parameters. 
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ice condition, that at every vertex there are exactly two arrows pointing in and two 
arrows pointing out. There are six possible configurations of arrows at each vertex, 
see the second row of pictures in Figure |2l hence the name of the model. 

To establish a correspondence between the three-colorings and the states of the 
six-vertex model, take a three-coloring, remove the boundary edges of the lattice and 
place arrows on the internal edges in accordance with the following rule. Consider 
four edges containing a fixed four-valent vertex of the lattice and four faces contain- 
ing these edges. The possible color combinations for such four-face sets are given by 
the first raw of the pictures in Figure |2l where r is one of the colors 0, 1, and 2. Visit 
the selected faces moving anticlockwise. If intersecting an edge we see that the color 
changes by +1 we place on the edge a pointing in arrow, if the color changes by —1 
we place a pointing out arrow. It is not difficult to get convinced that this rule is not 
contradictory, and we obtain the vertex configurations of the six-vertex model, see Fig- 
ure |2] The state of the six- vertex model, corresponding to the state of the three-coloring 
model given in Figured] is given in Figured It is also evident that the established cor- 



Figure 3: 



respondence is three-to-one, and if we fix the color of any face of the lattice it becomes 
one-to-one. 

In the present paper we show that similarity between the three-coloring model and 
the six-vertex model extends further. Namely, choose the weights of the six-vertex 
model as 



a(cp\fj) = a\(p\rj) 



sin(?//2 — (p) 



K<p\v) =P'(<p\n) 



sin tj 

701*7 ) = nr'(<p\y) = l > 



sin(///2 + (p) 
sin?/ 



(10) 
(11) 



see Figure |2] for the correspondence between the functions and the vertex configura- 
tions. The parameter r\ is called the crossing parameter. It is known that for r\ = 2n/3 
and the domain wall boundary conditions the partition function of the inhomoge- 
neous six- vertex model satisfies some simple functional equations HZHH, see Section |2] 
for details. We show that the partition function of the three-coloring model for appro- 
priate boundary conditions satisfies a similar equation. 

In the papers (ZHSl to obtain the functional equations for the six- vertex model the 
representation of the partition function via the Izergin-Korepin determinant was 
used. This representation is usually proved with the help of some recursion relations 
satisfied by the partition function which were found by Korepin IfTOH . In Section |2] we 
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obtain the functional equations directly from the recursion relations. In Section |3] we 
find recursion relations satisfied by the partition function of the three-coloring model 
for appropriate boundary conditions, which allow us to prove in this case functional 
equations similar to the functional equations of the six-vertex model. 

One can show that the statistical three-coloring model is a partial case of the cyclic 
SOS (solid-on-solid) model fTT| which in turn is a partial case of the eight-vertex SOS 
model |fl2|. Actually the weights of the eight- vertex SOS model depends on a param- 
eter 7/ analogous to the crossing parameter of the six-vertex model. The cyclic SOS 
model corresponds to the case when 7/ = 2n/L for some positive integer L, and the 
statistical three-coloring model corresponds to the case when rj = 2n/3. 

Recently it was shown by Rosengren that the partition function of the eight- vertex 
SOS model can be represented as a linear combination of determinants f!3|. In our 
case it is a linear combination of two determinants. One can use this representation to 
derive the functional equations, but we prefer to follow a direct way. 

The weights used in the paper ||T3l for t] = 2n/3 are different from the weights 
given by the formulas (H)-®. In the Appendix we show that these two sets of weights 
are connected by a change of the spectral parameter and a gauge transformation. 



2 Functional equations for six-vertex model from recur- 
sion relations 



In this section we consider the six-vertex model with the domain wall boundary con- 
ditions fllP], Here the boundary horizontal arrows point in and the boundary vertical 
arrows point out, see an example in Figure |U We label the horizontal lines by the 



— » — 1> 



-> — Q- 



4— «— X2 



Xi 



4— «— X3 



s> — *- 



X4 



#1 ^2 J/>3 ^4 



Figure 4: 



variables Xu * = 1, • • • , ft, and the vertical lines by the variables xpi, i = 1, . . . , n. The 
vertex at the intersection of the line labeled by Xi an d the line labeled by xpj acquires 
the spectral parameter Xi — tyj- The weight of a state of the model is the product of 
the weights of the vertices, and the partition function Z n ({x}; {</>}) is the sum of the 
weights of all possible states. 

The weights defined by the formulas (jlO|) and (flT) satisfy the Yang-Baxter equa- 
tions. Using this fact, one can show that the partition function is separately symmetric 
in the variables Xi and tyi, see, for example, ||T4llT5f . 
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Xp 1 Xp 2 $3 l/> 4 



Figure 5: 

Define the function 

n n n 

Fn({x}; W) = n sin te - xj) EI sin te - ty) n sin (^ - </v) z »({*}; (V 7 })- 

t,y=i ! /j = i /,y=i 

Our goal is to prove that for 7/ = 2zr/3 the function F n {{x}> {V 7 }) satisfies the func- 
tional equations 

2 

E F » ■ ■ ■ / Xfc + 2tts /3, . . . , = (12) 

s=0 

and the functional equations^] 

2 

Y^F n ({x};ip 1 ,...,ip k -2ns/3,...,ip n ) = 0. (13) 

s=0 

These equations can be used to solve some enumeration problems for the alternating- 
sign matrices. In particular, it is possible to reproduce [7] the refined enumeration of 
the alternating-sign matrices conjectured by Mills, Robbins and Rumsey ||T6l IT7| and 
proved by Zeilberger [fl8|. 

We start with the case when t] is arbitrary Consider the vertex on the intersection 
of the lines with the labels Xn and ip n . There are two possible configurations of this 
vertex, the second and fourth ones in Figure |2l They have the Boltzmann weights 
sin(j//2 — Xn +ipn)/ sin?/ and 1 respectively If we put Xn = tfn + ^/2, then only 
the lattice configurations corresponding to the second possibility will give a nonzero 
contribution to the partition function. Here the configuration of the vertices belonging 
to the bottom row and the left column become fixed, see Figure [5] for an example. 
The configurations of the remaining vertices correspond to the configurations of the 
(n — 1) x (n — 1) lattice with the domain wall boundary conditions. Collecting the 
weights of the fixed vertices, we come to the following recursion relation 



2 We write the equations fT3] | with a minus sign before 2ns/ 3 in the right hand side for similarity 
with the corresponding equations arising in the case of the three-coloring model. In the case of the 
six-vertex model we can use a plus sign as well. 



ZniXli- ■ -'Xn-l,Xn;ipl,- ■ • ,%-l f ^n)\ Xn= ^ n+T j/ 2 

n—1 n — 1 

= sin 2 " 2 " rj Y[ ^{Xi -xp n + rj/2)Y[ sin(i/?„ - ft + rj) 

i=l 2=1 

x Z„_i(^i, . . . ,x„-i;ipi, . . . ,fti-i). 

Actually, since the function Z n ({x}; {xp}) is symmetric in the variables Xi an d $i> we 
have n 2 recursion relations 



\x k =tpt+tl/2 



n 

.2-2w 



sin " rj Yl sm(xi - ft + rj/2) Y[ sin(ft^ - ft + 77) 



i=i i=i 



x Z„_i(ki,. . . .,#„; ft,. . .,ft,. . .,ip n ). (14) 

Now we start with the vertex on the intersection of the lines with the labels Xi an d 
ft,. In a similar way as above we obtain n 2 recursion relations 



Z n {Xl> ...,Xk,---,yn;ip\,---,lpl,---,lpn, 



n 

= sin 2-2 " t] Yl sm(xi - ft - rj/2) Y[ sin(ft -ipi-t]) 

z=l z=l 

i^k i£t 

X Zn-liXl, • • • ,Xk, • • • ,Xn,^\, • • • • • • ,%)■ (15) 

The recursion relations ((14)) and ((151) are valid for an arbitrary r\. Now assume that 
rj = 2k/ '3. Using the identity 

1 

sin(^) sin(^ + n/3) sin(tp + 2n/3) = - sin(3^) 



and the recursion relations ((T4|) and (flB) , we can prove that 



f n Ul/ • • • / Xh ■ ■ ■ / Xn) $\, ■ ■ ■ , Ipi, ■ ■ ■ , tp n ) \ Xk=t p e ± n /3 

= ± (_i)»-i 4 2 - 2m sin 3 - 2 "(27r/3) nsin[3(^ - Xi)\ flsin^ft - ft)] 

z=i z=i 

X ?n-\{X\> • • • >Xh • • • ,Xn)$\, ■ ■ ■ /tyi, • • • (16) 

Define the functions 

2 

s=0 

and prove by induction that S« ; fc({/t}; {ip}) = 0. Actually due to the skew-symmetry 
of the functions F n ({x}', I 1 / 7 }) ^ the variables it suffices to prove that 

s n ,i{{xh{f}) =o. 
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It is not difficult to see that 



2 

Si,iUi;V>i) = YL sin (^i - #1 + 27 */3) = o. 

r=0 



Assume now that S n _i r i({^}; {*/>}) = for some n > 1. The recursion relations (|T6)) 
give 

S n ,l(^l,- • •,^n-l/^«;'/ 7 l'- • • •' , / 7 «)l x „=^±7r/3 

= T (-l)"4 2 - 2 " sin 3 - 2M (27T/3) f[ sin[3(</> £ - X i)]fl si ^i ~ %)} 

i=l j=l 

X ^>n-\,\{X\> • • • >Xn-\)^\, • • • • • • ,^r?) = 0. 

It is clear that for the domain wall boundary condition any row or a column of 
a state has an odd number of 7-type and 7 / -type vertices. Therefore, we have, in 
particular, that 

Zn(Xi>--->Xn + n;{ip}) = {-l) n - 1 Z n {xi,---,Xn,{ip}). (17) 

With respect to the variable Xn the partition function Z n ({x}; {$}) is a trigonomet- 
ric polynomial of order less or equal to n — 1, and S n ,i({^}; {$}) is a trigonometric 
polynomial of order less or equal to 3n — 2. As we showed just above, the function 
S ni i({x}', {V 7 }) nas zeros at the points Xn — tf>£ =t n/3, t = 1, . . . , n. Moreover, by con- 
struction it also has zeros at the points Xn = Xl> ^ = 2,..., n — 1 and at the points 
Xn = tf>£/ t — l,...,n. The relation (|T7|) doubles all these zeros. Hence, we have 8n — 4 
zeros in the interval < Xn < 2n. This is possible only if S n/ i({x}', {V 7 }) = 0- Thus, 
we proved the equations ((T2j). The equations (|T3"1) can be proved in the same way. 



3 Functional equations for three-coloring model with do- 
main wall boundary conditions 

For the three-coloring model there is a natural analogue of the domain wall boundary 
conditions of the six-vertex model. Here if one starts with any boundary face and 
walks anticlockwise along the boundary then the color changes by +1 from face to 
face for the vertical boundaries, and by —1 for the horizontal boundaries. An example 
is given in Figured It is clear that in this case the number of face rows should coincide 
with the number of face columns. One can easily get convinced that the corresponding 
states of the six-vertex model satisfy the domain wall boundary conditions. 

We will consider the inhomogeneous case when the internal horizontal lines are 
labeled by the variables Xu z = 1, • • • ,n, and the internal vertical lines are labeled 
by the variables ipi, i = 1, . . . , n, see for example Figure [6l Here and henceforth n 
means the number of internal vertices in a row or in a column. With the vertex at the 
intersection of the line labeled by Xi and the line labeled by ipj we associate the spectral 
parameter Xi ~ We represent the total partition function Z n ({x}; {$}) as the sum 
of partial partition functions: 

z»({*};M)= Lz r n({xY,{f}), 

reZo 
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where r is the color of the left topmost face of the lattice. Note that the equalities 0-® 
imply that 

Zn + \{xh MIA, p) = Z" n ({xh + 27T/3, p). 

We denote the weight corresponding to the four-face configuration 



by W rs s O). As we noted above, these weights satisfy the Yang-Baxter equations 
which has the explicit form 

E w;7'(«p)w;;v)w t s ; v (<p -yi- tt/3) 

tez 3 

= EK>r(<p-?'-n/mY(<p')H s s '(cp) (18) 

tez 3 

and guarantees the integrability of the model. 

It is not difficult to get convinced that if weights Wjfg (<p) satisfy the Yang-Baxter 
equations ((181) , then for arbitrary non-zero constants C r gauge transformed weights 
W rs s 0), defined as 



rs C s , <P r {cp)& s ,(cp) rs ^ 

also satisfy the Yang-Baxter equations relation provided that 

<P r ((p) 



(19) 



&r{<p- <p' -7T/3) 



€v(<p')" 



The corresponding partition sums coincide in the case of the toroidal boundary condi- 
tions, and may differ by a factor in other cases. 
We apply the gauge transformation ((T9l> with 



C r = 1, <P r (q>) = C. 



l/12+<p/47T 



to the weights of the three-coloring model defined by the functions (OQ)-©. As the 
result we come to the weights defined by the functions 



U<P) 

lr{<p) 



0i(2n/3) ' 
1/2 9^71/3 + cp) 



# ( *° - 6 t (2n/3) 
e A (A + 2n(r + l/2)/3 + (p) 

9 A (X + 2nr/3) 
g 4 (A + 27r(r-l/2)/3- <p) 

6 4 (\ + 2nr/3) ' 



(20) 
(21) 
(22) 
(23) 



The corresponding partial partition functions for our boundary conditions are con- 
nected by the equality 



z£({*};W)=Il 







r+i 



-j(Xi - + Xn-i+1 ~ %-i+h 



®r+i(Xi ~ $i + Xn-i+1 ~ tfn-i+l) 



z;(k);W). 



Using a trick similar to that used for the six- vertex model ffT4l[T5| , one can show that 
the Yang-Baxter equations lead to the separate symmetry of Z r n ({x}; {ip}) in the vari- 
ables Xi an d i/ 7 /- 

Consider the vertex on the intersection of the line with the label Xn and the line 
with the label ip n . There are two possible four-face configurations for this vertex, the 
second and fifth ones in Figure |2l They have the Boltzmann weights &' r (Xn — tyn) and 
7r{Xn — tyn) respectively If we put Xn = $n + Tt/3, then only the configurations cor- 
responding to the second possibility will give a nonzero contribution to the partition 
function. Here the four-face configurations for the vertices belonging to the inter- 
nal bottom row and the internal left column become fixed, see Figure [7] for an exam- 
ple. The configurations of the remaining faces correspond to the configuration of the 
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Figure 7: 

(n — 1) x {n — 1) lattice with the domain wall boundary. Therefore, we have the fol- 
lowing recursion relation 



10 



Zn\Xl>- • • 'Xn-l/Xn'/flf ■ ■,%-!,%] 



Xn = 1pn + n/3 



= 7r(2n/3) n P'r+n-iiXi ~ %) Yl pr+n-i{fn ~ % + 71 / 3) 
i=l i=l 

x z r n _ 1 {xi,...,Xn-i;^\,- 

Using the definitions (|20)) - (j23|) and (0), we rewrite this relation as 



Z r n {Xl>- ■ • >Xn-l/Xn'/fh • ■ ■ ,tyn-\,tyn 

e A (X + 2n(r + n)/3) n ~ x 



ef- In (2n/3) 

n-l 



X 



e^lrVn-nk n *<* - * + ^ n * (*. - ft + M) 

X Z r n _ 1 (xi,...,Xn-i;il>l,-~,1pn-l)' 

Since the functions Z r n ({x}'r {V 7 }) are symmetric in the variables and we have 
actually n 2 recursion relations: 



e 2 " 2 "(27r/3) 



X 



jA+M^t-i^ n^' - * + ^ - ^. + ^/ 3 ) 

X Zr n _ 1 (xi,...,Xb...,Xn;ipl,---,$t,---,ll>n)- (24) 

Starting with the right topmost face and exploring again the symmetry of the func- 
tions Z^(x',ip) in the variables Xi and ty, we obtain another set of recursion relations: 



Zn (Xlr ■ ■ ■ > Xk, ■ ■ ■ , Xn) $1, ■ ■ ■ , ■ ■ ■ / $n) 



d\- lr - \2n/3) 



Xk=ipe-n/3 



X 



Finite - - 71/ 3) YlOltyt - % - ^71/3) 



i=\ 



i=i 



X Z r +\{x\,...,Xk,---,Xn)$\,---,$l,---,yn)- (25) 



The relations ((24)) and p5) are similar to the relations ((14)) and ((15)) of the preceding 



Section. Now we will find the analogue of the equations ((12)) and p3)). To this end we 
start with introducing the functions 



4 (A + 27r(r + n)/3) 

x n «iOc<-^) n hbi-ip)) n ^.-^(tow)- ( 26 ) 

i,/=l i,7=l 1,7=1 

It is not difficult to get convinced that 

e 1 (cp\p)e 1 {cp+ n/3\p)e l { ( p + 2n/3\p) = D(p)e 1 (3( P \p 3 ) / (27) 
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where D(p) is given by the relation 

D(P) = 38$^) • 

Using the equality ((27]) and the recursion relations ((24)), we see that 

Fn (Xl, ...,Xh---,Xn;ip\,---,yi,---,lpn\Kv)\ Xk =f e +n/3 

= (-l)"- 1 D 2w - 2 (p)^- 2 "(27r/3|p) n^(3(^ - ^)b 3 )fl^(3(^ - ^)lp 3 ) 

X fj|_i0tl/ • • • ,Xh ■ ■ ■ ,Xn) $1, ■ ■ ■ , h' ■ ■ ■ > A, p). (28) 
In a similar way the recursion relation ((25]) gives 

Kixi, •••,Xh---,Xn;fi,---,ip£,---,fn\h / p)\ Xk=xPe _ n/3 

= -(-l)"- 1 D 2 "- 2 (p)^- 2 "(27r/3|p) flOimi ~ Xi)\p 3 )fl0i(3(tpt - ^)\p 3 ) 

i=l /=1 

x F^lUi, • • -,Xh- • • ,Xn)^\, • ■ $1, • • ./</>n|A, p). (29) 
Now we define the functions 

2 

S n,k(Xl, ■ ■ ■ ,Xk, ■ ■ ■ >Xn) W) = Yu F n +S (Xl, ■ • • ,Xk + 2ns/3, . . . ,Xn) {</>})• 

s=0 

It appears that S r n k {{x}> I 1 / 7 }) = 0- To prove this fact we note first that it suffices to 
prove that S Q n 1 ({^}, {$}) = 0. All other cases can be reduced to this one using the 
skew-symmetry of the functions F£({x}/ {V 7 }) ^ n the- variables Xi an d the equality 

K + ~\{X), {<P}\K p) = mxh M|A + 27T/3, p), 
which follows from the relation 

^ S MX\^\ X 'V) = fi#(?|A + 27r/3,p). 
We prove the equality S Q n i{{x\r {$}) = by induction on n. It is easy to see that 

fr( , \ 0i(Xi ~ tpi)9j(X + Xi~ipi + 27r(r + 1/2) /3) 
tiKXl'WD 9 i (X + 2n(r + l)/3)9 A (X + 2nr/3) 

Hence, the function S° ^{xi'i Vi) considered as a function on cp = xi — ^i has the form 

0l(<P)04(A+<p+7T/3) 



SIM 



6 4 (X + 2tt/3)6 4 (X) 
e 1 ((p + 2n/3)9 4 (X + ( p + 5n/3) 6 x {(p + 4tt/3)6> 4 (A + <p + 3n) 
4 (A + 47r/3)0 4 (A + 27r/3) 9 A (X + 2n)e A (X + 4zr/3) ' { ' 
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All terms of the sum in the right hand side of the equality ((30|) have the same period- 
icity factor associated with the periods n and tit. Dividing this sum by its first term, 
we obtain a doubly periodic function. It is not difficult to verify that this function has 
no poles, hence, it is a constant At cp = —n/3 the function S^Jcp) turns to zero, but 
the first term of the sum in the right hand side of the equality (|30|) is not zero at this 
point. Thus, the constant in question is zero, and S® 1 (xi'/tf>i) = 0. As we noted above, 
this means that S\ 1 (xi ; =0 for any r. 

Now we assume that S r * \{{x}> {$}) — for some n > 1. It follows from the 
recursion relations (|28l) that 



S n,l {Xl>--->Xn-l> Xn> $\, • • • ■ ■ ■ ,1pn\K p) 



n-1 



4r- 1 D 2 "- 2 (p)^ 2 "(27r/3|p)n^(3(^-^OIp 3 )n i( 3 (^-^lr 



i=l i=l 

x S r n _ xl {xi, ■ ■ ■ ,Xn-Y,ipi, ■ ■ ■ ■ ■ ■ ,$n\Kv) = °' 



and the recursion relations (|29|) give 

Sn i (XV Xn-l, Xn) Ql, ■ ■ ■ , $1, • • • , $n\K P) 



Xn=ip£-n/3 



= -(-l)"- 1 D 2 "- 2 (p)^- 2 "(27T/3|p) n»l(3(^ " ^)IP 3 ) fl^(3(^-^)lp 3 ) 

i=l i=l 
X 5 n + -l,l{Xl> ■ ■ ■ >Xn-l^l, ■ ■ ■ ,^Pi, ■ ■ ■ ,^n\Kp) = 0. 

Recall that for any state of the six-vertex model with the domain wall boundary 
conditions the number of 7-type and 7 / -type vertices in any column and in any row 
of vertices is odd. It follows from the correspondence between the states of the three- 
coloring model and six-vertex model, described in Section [TJ that for any state of the 
three-coloring model with the domain wall boundary conditions the total number of 
7-type and 7'-type vertices in any column and in any row of the internal vertices is 
odd as well. Having in mind the equalities 

0i((p ± n\p) = -6i(<p\p), 4 (<p ± n\p) = 9 4 ((p\p), 

we conclude that any summand of the partition function Z r n ({x}', I 1 / 7 }) have the same 
periodicity factor associated with the period n with respect to any variable Xi- The 
same is true for the summands of the functions F„({x}; an d S r n k {{x}> I 1 / 7 })- 
Using the equalities 

di(q> + ttt\p) = -p- l e-np- 2icp 6i(cp\p), 6±{<p + nr\p) = -p" 1 exp -2 ^ e A {cp\p), 
we see that 

Wf/icp + TZt) = -p-V 2l ^wr/(<p), 

where 

Y = e i(f-l)[7r(r+l)/3+A]r 
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Now it is possible to show that any summand of the partition function Z r n ({x}; {ip}) 
have the same periodicity factor associated with the period tit with respect to any 
variable Xi- The same is true for the summands of the functions F r n {{x}> {$}) an d 

sk({*};M). 

Any summand of the function S r n {V 7 }) contains a product of 3n — 1 theta 

functions whose arguments contain the variable Xn- Since all summands have the 
same periodicity factors, if we divide S r nl ({x}', {V 7 }) by any of its summands we will 
obtain a doubly periodic function of Xn- Every elliptic theta function has one and only 
one zero in a cell with the corners £, £ + n, £ + n + tit, £ + T. Therefore, the order of the 
function obtained after division of the function S r n x ({#}; {$}) by any of its summands, 
is less or equal to 3n — 1. From the other hand, as we showed above the function 
S r nl ({x}; {V 7 }) is zero ror PCw = ^ =t 27T/3, £ = 1, . . . , n. Moreover, it follows from 
the definition of the functions S r n k ({x}; {</>}) that S r n l ({x}; {$}) is zero for Xn = Kb 
i = 2, . . . , ft — 1 and for # n = tpi, £ = 1, . . . , n. Hence, the function Sj^ {i/ 7 }) has 

at least 4n — 2 zeros, and since 4n — 2 > 3n — 1, we have {i/?}) = 0. 

Thus, we proved that 

2 

E K +S {Xl, ■■■,Xk + 2tts/3, . . . ,Xn) W) = 0, 

s=0 

where the function V r n ({x}> I 1 / 7 }) is related to the partition function of the model via 
the equality (|26j) . In the same way we can show that 

2 

E F n +S ({XY> lpl,...,lpk- 27TS/3, . . . , fa) = 0. 

s=0 
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Appendix 

After the substitutions A — > A + tit/ 2 and cp — > —cp— n/3, using the relation 

9 4 ((p\p) = i/^e-^fliO - TtT/2\p), 
we transform the weights ©-© to the weights 

<f)=<(f)=^' t " H ^^ ) , (3D 

M?) = IH(9) = -C' 2+MllK ^3), (32) 

lrM = e'f 2^ 9.(A + 2^r/3-y) 

^ rf/2n fl 1 (A + 27rr/3) ' 1 ' 
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V (cp) = e ^f ^ ft(A + 27rr/3 + y) 

e 1 (A + 2tit/3) ' 1 j 

b r-l 

where 

01 (A + 2n(r - 1)/3)9 1 (A + 2zr(r + l)/3) 
02(A + 2m73) 

These weights satisfy the Yang-Baxter equations of the form 

E w#>)w;;V)< s '(<p-<?>') = E Wr>Hv-v')</W)Hi{<p) 

tez 3 tez 3 

and the gauge transformations are again given by the formula (fl9~1) , where the func- 
tions <f> r (<p) satisfy now the relation 



More explicitly the gauge transformations (fl9[) are described by the formulas 

= — &K<pf — = — ^K<pf — 

P ^ C ,-T * ^(f )<2> r+I (<p) PrWJ ' PrW C r+T *^)* r+I (^) PAV> ' 

Applying the gauge transformation (fl"9)) with 

Cr = e i7T/ V /2 (A + 27rr/3), 4> r (p) = ^<p/2^f/^ 



to the weights (|3T|) - (|34|) we obtain 

_ , 0i(2tt/3 + <p) , 0i(27r/3 + (p) 

^ = ^3) - ^ = 0i(2zr/3) ' 

^ (y) _ g i( A + 2 ^~l)/3) 0i(y) 0, (9) _0i(A + 27r(r + l)/3) 0i(?) 



0i(A + 2m73) 01(271/3)' rrVr; 0i(A + 2m73) 0i(2tt/3)' 
. , 0i(A + 27rr/3- <p) . X (A + 2zrr/3 + <p) 

= 0i(A + 2zrr/3) - ^ = 0i(A + 2zrr/3) ' 

These weights coincide with the weights used in the paper ffl3|. 
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